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Abstract. We review our studies of elastic lattices driven by an external force F in the presence of 
random disorder, which correspond to the case of vortices in superconducting thin films driven by 
external currents. Above a critical force Fc we find two dynamical phase transitions at Fp and Ft, 
with Fc < Fp < Ft. At Fp there is a transition from plastic flow to smectic flow where the noise is 
isotropic and there is a peak in the differential resistance. At Ft there is a sharp transition to a frozen 
transverse solid where both the transverse noise and the diffussion fall down abruptly and therefore 
the vortex motion is localized in the transverse direction. From a generalized fluctuation-dissipation 
relation we calculate an effective transverse temperature in the fluid moving phases. We find that the 
effective temperature decreases with increasing driving force and becomes equal to the equilibrium 
melting temperature when the dynamic transverse freezing occurs. 

INTRODUCTION 

The study of the collective motion of vortex lattices in superconductors has brought 
new concepts in the non-equilibrium statistical physics of driven disordered media 
[1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13]. The prediction by Koshelev and Vinokur [1] 
of a dynamical phase transition upon increasing drive, from a fluid-like flow regime 
[2, 3, 4] to a coherently moving solid [1], has motivated an outburst of theoretical 
[5, 6, 7, 8], experimental [9, 10, 1 1], and simulation [12, 13] work. The relevant physics 
of the high velocity driven phase is controlled by the transverse displacements (in the 
direction perpendicular to the driving force) [5], leading to a new class of driven systems 
characterized by anisotropic spatial structures with transverse periodicity [5, 6, 7]. 
These moving anisotropic vortex structures have been observed experimentally by Pardo 
et al. [11], and their different features have been studied in simulations [12, 13]. 

Koshelev and Vinokur [1] defined a "shaking" temperature from the fluctuating 
force felt by a vortex configuration moving in a random pinning potential landscape. 
This lead to their prediction of a dynamic phase transition when T^h equals the equilib- 
rium melting temperature of the vortex system [1]. This dynamic transition separates the 
liquid-like phase of vortices moving at low driving forces from the "crystalline" vortex 
lattice moving at large forces. However, later work [5, 6, 7] has shown that the perturba- 
tion theory used in [1] breaks down and that the vortex phase at large velocities can be an 
anisotropic transverse glass instead of a crystal. In spite of this, the shaking temperature 
introduced in [1] has been a useful qualitative concept, at least phenomenologically. 

Cugliandolo, Kurchan and Peliti [14] have introduced the notion of time-scale de- 
pendent "effective temperatures" T^ff from a modification of the fluctuation-dissipation 
theorem (FDT) in slowly evolving out of equilibrium systems. Teff is defined from the 



Nonequilibrium structures and dynamic transitions in driven vortex lattices with disorderFebruary 6, 2008 1| 



slope of the parametric plot of the integrated response against the correlation function of 
a given pair of observables when the latter is bounded or, equivalently, it is the inverse 
of twice the slope of the parametric plot of the integrated response against the displace- 
ment when the correlation is unbounded. This definition yields a bona fide temperature 
in the thermodynamic sense since it can be measured with a thermometer, it controls the 
direction of heat flow for a given time scale and it satisfies a zero-th law within each time 
scale. Teff was found analytically in mean-field glassy models [15] and it was success- 
fully studied in structural and spin glasses, both numerically [16] and experimentally 
[17], in granular matter [18] and in weakly driven sheared fluids [19]. 

In this article, we will review our work on the dynamical regimes and nonequilibrium 
transitions of the driven vortex lattices as a function of the driving force [13], and show 
how an adequate definition of an effective temperature (based on Ref. [14]) can help to 
understand these phenomena [20]. 



IMODEL 

The equation of motion of a vortex in position T; is: 

= - E V,C/,(r,-,-) V,-C/^(r,p) + F+ Ci(t) (1) 

where rij = \ri — rj\ is the distance between vortices rip = |r, — r^l is the dis- 
tance between the vortex / and a pinning site at r^, t] = ^'^2p^^ is the Bardeen- 

Stephen friction and F = x z is the driving force due to an applied current J. 
In two-dimensional superconductor the vortex-vortex interaction potential is logarith- 
mic: Uv{r) = — Avln(r/A), with Ay = 4>Q/8;rA [3]. The vortices interact with a random 

uniform distribution of attractive pinning centers with Up{r) = — A^e^^''/'''')^ with rp 
being the coherence length. We normalize length scales by rp, energy scales by Ay, 
and time is normalized by T = rir^/Ay. The effect of a thermal bath at temperature 

T is given by the stochastic force Cii^), satisfying {<^^{t)) = and {C,^{t)C,^ {t')) — 
lTi]kBT5{t — t')5ij5^^i. We consider A^^ vortices and Np pinning centers in a rectangular 
box of size x Ly, and the normalized vortex density is = Nyr^/LxLy = Br^/<^Q. 
Moving vortices induce a total electric field E = x z, with v = ■^'Li'^i- We use pe- 
riodic boundary conditions and the periodic long-range logarithmic interaction is eval- 
uated with an exact and fast converging sum [21]. Typical simulation parameters cor- 
respond to vortex densities of = 0.05 — 0.12 in a box with Lx/Ly = \/3/2, with 
A^v = 64 — 784, pinning strengths of Ap/A^ = 0. 1 — 0.35, and densities of pinning cen- 
ters Hp > Hy (dense pinning Hp > Uy, is typically realized in experimental samples). The 
equations are integrated with a time step ofAf = 0.01— O.It and averages are evaluated 
in 30000 - 80000 integration steps after 2000 - 30000 iterations for equilibration (when 
the total energy reaches a stationary value). 
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FIGURE 1. Examples of steady states for different dynamical regimes. Upper panel: Vortex trajectories. 
Middle pannel: Average structure factor 5'(k). Lower pannel: Voronoi construction for each given vortex 
structure. 



DYNAMICAL REGIMES 

We start by showing the vortex trajectories and their translational order in the different 
steady state regimes. In the upper panel of Fig. 1 we show the vortex trajectories {r, (?)} 
for typical values of F by plotting all the positions of the vortices for all the time iteration 
steps. We also study the time-averaged structure factor 5(k) = (|^X^,exp[«k ■ r,(?)] p), 
which is shown in the middle pannel of Fig. 1. The number of deffects in the lattice 
structure can be obtained from a Voronoi construction. This in shown in the lower pannel 
of Fig. 1, where the lattice defects are shown in gray. In Fig. 2(a) we plot the average 
vortex velocity V = {Vy{t)) — i^'Li^)^ in the direction of the force as a function of 
F and its corresponding derivative dV /dF. Below a critical force Fc all the vortices 
are pinned and there is no motion, V — 0. Above Fc, we can distinguish three different 
dynamical regimes: 

( i) Plastic flow: Fc < F < Fp. At Fc vortices start to move in a few filamentary channels 
[3]. A typical situation is shown in Fig. 1, where a fraction of the vortices are moving in 
an intricate network of channels. As the force is increased a higher fraction of vortices 
is moving. In this regime, vortices can move in the transverse direction (perpendicular 
to F) through the tortuous structure of channels [4]. We see that the corresponding 5(k) 
only has the central peak showing the absence of ordering in this plastic flow regime 
[2, 3, 4] and that there are a large number of defects. 

(ii) Smectic flow. Fp < F < Ff. We find a peak in the differential resistance, dV /dF, 
at a characteristic force Fp. At F = Fp we see that all the vortices are moving in a seem- 
ingly isotropic channel network with maximum interconnectivity. In the experiment of 
Hellerqvist et al [10] the value of Fp was taken as an indication of a dynamical phase 
transition. In fact, we find that above Fp a new dynamic regime sets in. In this case, as 
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we show in Fig. 1, all the vortices are moving in trajectories that are mostly parallel to 
the force, forming "elastic channels". Two small Bragg peaks appear in the structure 
factor along the Ky = axis, which correspond to Gi = (±2;r/ao, 0)- This is consistent 
with the onset of "smectic" ordering [6] in the transverse direction with elastic channels 

— 1/2 

separated by a distance ~ ao = «v -In this regime the transverse motion consists in 
vortex jumps from one channel to another, resembling "thermally" activated transitions 
induced by local chaos. The rate of these transitions decreases with increasing force. In 
Fig. 2(b) we plot the magnitude of the Bragg peaks at Gi, Ss = S{G\), corresponding to 
smectic ordering (Ky = 0), and the other neighboring peaks at G2 = ±2;r /ao ( 1 /2, V3/2) 
and G3 = ±2n/ao{ — 1/2,^/3/2) , Si = {S{G2) +S(G3))/2), corresponding to longitu- 
dinal ordering (Ky 7^ 0). We see that above Fp the intensity of the smectic peak Ss starts 
to grow and Ss ^ Si, while below Fp the spatial structure is isotropic, Ss — Si <^ 1. The 
Bragg peak heights depend with system size as 5(G) ~ A^v where <7g = means 
long-range order (LRO), < ctg < 1 rneans quasi long-range order (QLRO) and (7g > 1 
means short-range order (SRO). We have found that ctg > 1 in this regime: there is only 
short range smectic order [13], and thus this phase corresponds to a fluid. In this sense 
the transition at Fp is a dynamic transition in the flow. 

Frozen transverse solid: F > Ff. At a new characteristic force Ft, the jumps 
between channels suddenly stop and vortex motion becomes frozen in the direction 
perpendicular to F. An example for F > is shown in Fig. 1 where we see well 
defined elastic channels parallel to F. We see that new peaks appear in 5(k) in directions 
with Ky 7^ 0, like G2, G3, showing that there is some longitudinal ordering between 
the channels. These extra peaks are smaller than the smectic peaks, and 5(k) is very 
anisotropic. We note in Fig. 2(a) that Ft is the point where the noisy behavior in 
dV /dF ceases. A similar criterion was used by Bhattacharya and Higgins to define 
their dynamical phase diagram [9]. In Fig. 2(b) we see that in Ft there is an increase 
in the longitudinal ordering Si, and both Ss and 5/ tend to saturate at an almost constant 
value for F ^ Ft. Furthermore, we have we found that there is QLRO with a value 
of 0.5 < aG < 0.7 [13]. We also see that in this phase there are very few defects, 
corresponding to dislocations oriented in the direction of the force. 

DYNAIVIICAL TRANSITIONS 

An understanding of the dynamical transitions can be obtained from studying the tempo- 
ral behavior of the system. It has been observed experimentally that the longitudinal volt- 
age can show low frequency noise. This voltage noise reaches a very large value above 
the critical current, which has been attributed to plastic flow, and then the noise decreases 
for large current. In addition, even when the total dc transverse voltage (V^) = ( ^ ^) 
is zero, it can also have fluctuations and noise. In fact, it is easy to understand that this 
transverse noise will be closely related to the wandering and wiggling of the plastic 
flow channels and to the jumps between elastic channels in the smectic phase. We have 
calculated in [13] the power spectrum of both the longitudinal voltage and the trans- 
verse voltage. We have found [13] that near the critical force, the longitudinal noise, 
Py, is larger than the transverse noise, Px- At Fp the voltage noise becomes isotropic. 
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FIGURE 2. (a) Velocity-force curve (voltage-current characteristics), + symbols. dV/dF curve (dif- 
ferential resistance), continuous line, (b) Intensity of the Bragg peaks. For smectic ordering, Ky = : 
triangles; for longitudinal ordering, Ky ^ 0, asterisks, (c) Diffusion coefficient for transversal motion, D^. 



Px — Py. This is the point where we have seen the highest interconnection in the chan- 
nel network. The coincidence of isotropic noise, onset of short-range smectic order and 
peak in differential resistance suggest that there is a dynamic transition at Fp. Above 
Fp, the onset of elastic channels and smectic ordering reduces the longitudinal noise, 
while the transverse noise remains large due to the "activated" jumps between elastic 
channels. At Ft the transverse noise falls abruptly, by several orders of magnitude. This 
corresponds to a freezing transition of vortex motion in the transverse direction. Above 
Ft there are no more vortex jumps between elastic channels. The low frequency noise 
can be closely related to diffusive motion for large times. We have analyzed the average 
quadratic displacements of vortices in the transverse direction from their center of mass 
position {Xcm{t),Ycm{t)) as a function of time. We define Wx{t) = ^^Y.i[^i{t) —Xi{0)]^, 
where Xi{t) — Xi{t) —Xcm{t). We have found that the vortex motion in the fluid regimes, 
for Fc < F < Ft, is diffusive in the transverse direction, with Wx{t) ~ Dxt. In Fig. 2(c) 
we show the behavior of the transverse diffusion coefficient D^. The transverse diffusion 
is maximum at the peak in the differential resistance, Fp, and it has a clear and abrupt 
jump to zero at Ft indicating the transverse freezing transition. It is interesting to note 
that melting transitions also show a jump in the diffusion coefficient. 
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EFFECTIVE TEMPERATURE 



To study the fluctuation-dissipation relation (FDR), we proceed in a similar way as in 
simulations of structural glasses [16] and consider the observables: 

Mf) = wL '4 (0 ; 5m (0 - 1 ^4 (0 , (2) 

where 5/ = — 1, 1 are random numbers with sj = and Jjs] = 5ij, and = — Rem 
with ^ = x,y and Rem being the center of mass coordinate. Taking F = Fy we study 
separately the FDR in the transverse and parallel directions with respect to F. The time 
correlation function between the observables A^i and is 

C^{tM=Wt)BM} = ^I.irf{t)rf{tQ)) , (3) 

smce the rf are independent of the Si. The integrated response function Xpi for the 
observable is obtained by applying a perturbative force ^ = esifi (where fi =x,y) 
at time to and keeping it constant for all subsequent times on each vortex: 



1 



X^tit,to) = lim - {A^{t))^ - (A^(0)e=o 



(4) 



We then define a function, T^^^(t,to), by the relation, 

Xnit,to) = ——^——A^,{t,to) , (5) 

where A^(f,ro) = w.^fZiilrf (t) - r';' {to)\^) = C^{t,t) + C^i{to,to) - 2C^i{t,to) is the 
quadratic mean displacement in the direction of fi. For a system in equilibrium at 
temperature T the FDT requires that T^^^ = T^^^ = T. In a nonequilibrium system, like 
the driven vortex lattice with pinning, the FDT does not apply. Since we are interested 
in the stationary states reached by the driven vortex lattice, where aging effects are 
stopped [19], then all observables depend only on the difference t — to, if we choose 
to long enough to ensure stationarity. From the parametric plot of vs. A^(?) one 
defines the effective temperature T^^^ (t) using Eq. (5), provided T^^^ (t) is a constant in 
each time-scale [14]. 

We have studied [20] the transverse and longitudinal FDR for the moving vortex 
lattice as a function of driving force, F, solving the dynamic equations for different 
values of Ap, n^,, and T within the fluid regimes, Fc < F < Ft. We have obtained [20] 
for each F parametric plots of the transverse quadratic mean displacements, Ax{t) as a 
function of the integrated transverse response, Xx{t)- We found that the equilibrium FDT 

P 



does not apply in general but two approximate linear relations exist for Ax{t) < O.OSr^ 



r2 
P' 

short displacements region corresponds to the bath temperature, T, meaning that the 



and for Ax{t) > Vp, with a non-linear crossover between them. We have found that the 
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FIGURE 3. Transverse effective temperature T^^^ vs voltage V for Ap = 0.35, T = and tiy = 0.07, 
using the diffusion relation (filled diamonds) and a generalized Kubo formula (open diamonds). Inset: 
Transverse diffusion constant DxysV. Dashed lines indicate the transverse freezing transition at F =Fi 
and the dash-dotted Une indicate the melting temperature of the unpinned system « 0.007. 

equilibrium FDT applies in the transverse direction only for short times, t <^ rp/v. For 
the large displacements region we have obtained a different slope corresponding to an 
effective transverse temperature T^ff{T) > T. Furthermore, when comparing the results 
for different T , we have found T^^^{T) ^ T^^^{0) + T. On the other hand, a similar 
analysis of the FDR for the longitudinal direction does not show a constant slope for 
large displacements or long times such that t > rp/v. 

In Fig. 3 we show the transverse effective temperature obtained in Ref. [20], T^ff, 
for r = as a function of voltage (i.e., average velocity, V). We observe that T^ff is 
a decreasing function of V. The most remarkable result is that T^^^ reaches a value 
close to the equilibrium melting temperature of the unpinned system, ~ 0.007 [22], 
when the system approaches the transverse freezing transition at F = Ft (obtained from 
the vanishing of the transverse diffusion D^, shown in the inset). We have found that 
Tgjf Tm when F ^ Fj for different values of pinning amplitude, Ap, and vortex density, 
riy, even when Ft depends on Ap, 

CONCLUSIONS 

In conclusion, we have obtained evidence of two dynamical phase transitions. The first 
transition at Fp is the point of isotropic noise and maximum transverse diffusion and 
corresponds to the observed peak in the differential resistance. The second transition 
at Ft is a freezing transition in the transverse direction, where the transverse diffusion 
vanishes abruptly. We have been able to define an effective temperature in the moving 
fluid phase, using the thermodynamically adequate definition of [14]. The fact that we 
find that dynamic freezing occurs when T^ff{T) — T^, clearly indicates that there is a 
nonequilibrium phase transition at Ft. 
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